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We use ultra-high-speed video imaging to look at the initial contact of a drop impacting 
onto a liquid layer. We observe experimentally the vortex street and the bubble-ring 
entrapments predicted numerically, for high impact velocities, by Thoraval et al [Phys. 
Rev. Lett. 108, 264506 (2012)]. These dynamics occur mostly within 50 /is after the first 
contact, requiring imaging at 1 million frames/sec. For a water drop impacting onto a thin 
layer of water, the entrapment of isolated bubbles starts through azimuthal instability, 
which forms at low impact velocities, in the neck connecting the drop and pool. For 
Re above about 12 000, up to 10 partial bubble-rings have been observed at the base 
of the ejecta, starting when the contact is ~ 20% of the drop size. More regular bubble 
rings are observed for a pool of ethanol or methanol. The video imaging shows rotation 
around some of these air cylinders, which can temporarily delay their breakup into micro- 
bubbles. The different refractive index in the pool liquid reveals the destabilization of the 
vortices and the formation of streamwise vortices and intricate vortex tangles. Fine-scale 
axisymmetry is thereby destroyed. We show also that the shape of the drop has a strong 
influence on these dynamics. 



1. Introduction 

The impact of a drop onto a pool surface has been studied for over a century, but 



2003) have 



revolutionary improvements in high-speed video technology (Etoh et al. 
recently opened up this canonical geometry to renewed scrutiny. This applies especially 
to the earliest contact between the drop and pool, where intricate details have emerged 
and play a crucial role during air entrapment and splashing ( Yarin [2006 Thoroddsen et 
~ n[2QQ8l ). 



The impact of a drop always entraps a bubble under the center of the drop, as a disk 
of air is produced by the lubrication pressure and rapidly contracts into a bubble at 



the center ( 


Thoroddsen et al. ||2003| Liow & Cole ||2007 


Korobkin et al. 


2008 


Mani et 


al. 1120101 |H 


icks & Purvis |2010t Driscoll & Nagel |2011 


|Kolinski et al. 


2012 


|van der 


Veen et al. 


2012). Following this central air disk entrapment on a liquid pool, the outer 



contact forms a neck, which emits an ejecta sheet for sufficiently large Reynolds numbers 



( Thoroddsen 


2002 Weiss 


& Yarin 


1999 Davidson 


2002 Josserand & Zaleski 


2003 


Howison et al. 


2005 ). These ejecta are the source of the finest spray droplets (Thoroddsen 


et al. 2011 Zhang et al. 


2012), which is of relevance to numerous processes, such as 



combustion and aerosol formation. 

However, at even larger impact energy, these ejecta give way to random splashing of 
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Figure 1. (a) Imaging setup. We use backlight imaging of the drop impact from below, through a 
glass plate. The drop acts as a lens focusing the illumination to a limited area on the observation 
window. Several lights were used in some cases to obtain a larger illu mina ted area, (b) Camera 
viewing area shown in bottom view, corresponding to the area in Fig. 10 'c). 



small droplets, see Thoroddsen (2002). Numerical simulations by Thoraval et al. (2012) 



have shown that the base of the ejecta can become unstable, bending up and down as the 
free surface sheds alternate sign vortex rings into the liquid and often entrapping bubble 
rings. These bubble rings alternate between the top and bottom sides of the ejecta. This 
regime is the focus of the current investigation. 



Very recent experiments (Castrejon-Pita et al. 2012) have used side-view and laser 



induced fluorescence to verify the presence of the von Karman street for conditions similar 



to those in Thoraval et al. (2012). Herein we show the flrst experimental observations 



of the formation of the bubble tori. 

The main axisymmetric features of the vortex street and bubble rings entrapments are 
observed experimentally. However, three-dimensional effects rapidly break the symmetry. 
Herein, we show that even at rather modest impact velocities, azimuthal instabilities can 
appear in the neck between the drop and the pool. Imaging using two different liquids 
also reveals the shedding of streamwise vortices and their intricate dynamics, similar 
to three-dimensional instabilities of the cylinder wake ( Williamson 111996 ), or the shear 



layer (Lasheras & Choi 1988). These intricate structures have perhaps escaped earlier 



experimental notice as they develop in a sub-millimeter region and evolve in less than 
50 /is. 



2. Experimental Setup and Numerics 

2.1. High-Speed Video Imaging 

In this work we image drop impacts onto shallow pools through a bottom glass plate 
(Fig. [T]). Limited imaging (only in Fig. 16(a,b)) was done from the side above the pool 
surface. We use identical water drops in the entire study, while changing the composition 
of the pool liquid. The pool liquids tested are water, ethanol and methanol which are 
all highly miscible with the water drop. The liquid properties are given in Table [1] The 
difference in refractive index between the water drop and the ethanol or methanol pools 
allows us to image the flow structures as they distort the interface between the two 
liquids. 

The use of shallow pools, or thin fllms, is dictated by the need to change the pool liquid 
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Liquid 




M [cP] 


ly [cSt] 


a [dyne/cm] 


n 


Distilled water 


0.996 


1.004 


1.008 


72.1 


1.333 


Ethanol 


0.789 


1.19 


1.51 


23.2 


1.363 


Methanol 


0.793 


0.593 


0.748 


22.5 


1.339 



Table 1. Properties of the different liquids used in the pool. Here p is the liquid density; /i is 
the dynamics viscosity; u is the kinematic viscosity, n the refractive index at A = 532 nm, and 
a the surface tension. The drop is always water. 



following every impact as well as by the optical setup, where the limited focal distance 
of the long-distance microscope rules out bottom views through deep pools. 

The pool depths 5 were varied from about 25 /im to 1 mm. Here we use a long- 
distance microscope for magnifications up to about 15 for maximum pixel resolution 
of ~ 4.1 /im/px, when using the Shimadzu Hypervision CCD video camera (Etoh ei| 
\aL ||2003 ), at frame rates up to 1 million fps. Some of the imaging was also done at 



a lower frame rate with a Photron SA5 CMOS camera, with a magnification up to 
10 and maximum pixel resolution ~ 2 /im/px. Using thin bottom layers restricts the 
vertical motion of the interface between the drop and the pool liquid during the impact, 
thereby making well- focused imaging easier with the limited focal distance. For further 
optical/triggering details see Thoroddsen et al (2012|). 



The drop is pinched from a 3 mm nozzle, to produce an effective drop diameter of D = 
{DyDf^y^^ = 4.67 mm, where Dy and are the instantaneous vertical and horizontal 
diameters. We characterize the impact conditions by the Reynolds number Re, the Weber 
number We and the splashing parameter defined as: 

pDV pDV'^ , — 

p a 

where p, p and a are respectively the density, dynamic viscosity and surface tension of 
the drop liquid, and V the drop impact velocity. 

The drop velocity V was characterized in a separate set of experiments. It was then 
modeled by the velocity of a sphere experiencing constant drag: 

V = VT^/l - exp {-2g {h-D- ho) /VS), (2.1) 

where gravity is g = 9.81 m/s^, with the fitting parameters Vt = 9.11 m/s and ho = 
2.1 mm. Vt corresponds to the terminal velocity of the drop, and ho to the pinch-off 
length of the drop when it separates from the nozzle. Here h is the measured distance 
from the nozzle tip to the undisturbed pool surface, whereas the adjusted height is defined 
diS H = h - D - ho. Figure |2ja) shows that the measured values of V are less than 0.8% 
away from the formula, for our impact heights 2.5 cm < H < 55 cm. This estimate of Vt is 
slightly higher than the experimental observations of |Gunn fc Kinzer] ( |1949 i|, which could 
be due to the drop oscillations before reaching a final oblate shape. We can calculate the 
falling time of the drop from the falhng height as: t = (VT/^)argch [exp (^gH/Vr^)]. 

As D is larger than the capillary length for water, Ic = \/o'/{pg) = 2.7 mm, the water 
drop shows large oscillations that can affect the details of the impact dynamics (see 
Fig. [2|b)). The axisymmetric vertical oscillations of the drop can be estimated by the 
dominant mode ( Rayleigh] ( |1879| ) , Lamb (1975} § 275)), giving a radius: 



R{t, 0) = Ro [1 + acos (uot + (j)) P2 (cos (9)] , (2.2) 
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Figure 2. (a) Drop velocity V vs. falling height H . The green line corresponds to \/2gH, while 
the red line is our fitting equation, (b) Typical drop shapes in air: prolate {a> 1, H = 31.2 cm, 
Re = 11 300), oblate (a < 1, if = 41.8 cm, Re = 12 900) and close to spherical {a c:^ 1, H = 44.5 cm, 
Re = 13 300). The scale bar is 2 mm long. 



where P2 (^) = (3x^ - l) /2 is the Legendre polynomial of degree 2, and is the polar angle 
in the spherical coordinate system. The aspect ratio between the vertical and horizontal 
diameters of the drop can thus be written: 

Dy 1 + a cos (cjt + ^) 
Dh 1 - f cos {ujt + (j)) 



(2.3) 



We determine a, u and (j) as fitting parameters: a = 0.162, / = a;/(27r) = 33.2 Hz and 
(f) = -132°. The oscillation frequency is only 2% lower than the inviscid theoretical value 
Id = (4/7r)>/cr/(pI)^) = 33.9 Hz. The typical time of bubble-ring entrapment, 50 /is, 
is only 0.17% of the oscillation period. Therefore the drop shape can be considered 
frozen during the entrapment. This fitting is then used to get the aspect ratio from the 
falling height in the experiments. We have neglected viscous damping of the dominant 
mode in this estimate of drop oscillatio ns. The chara cteristic time of this damping can 
be estimated as r = D'^/{20i') 
studied here 0.35 s). 



1.08 s (Lamb 



1975 



355), In the overall falling time 



viscous effects can be estimated to reduce the amplitude of the 
dominant mode by 27%. It is therefore too short to damp the oscillations significantly, 
as is observed in Fig. [3j 

2.2. Numerical method 



We use the open-source software Gerris ( |http: / / gfs.sf.net Popinet (2003, 2009); Ag- 



baglah et aFI ( |2QllD ), using the Volume- Of- Fluid method, to perform axisymmetric sim- 
ulations of the drop impacts. The liquid from the drop and the pool are identified with 
different markers (drop: red, pool: blue, air: light green). The adaptive mesh is refined 
dynamically based on the distance to the interface, vorticity magnitude and geometric 
conditions. The interface is refined uniformly at the maximum level in the simulations. 
The bubbles and droplets with area less than 10 cells are removed during the computa- 
tion, as their dynamics cannot be captured accurately. It represents an effective cutoff 
diameter of Dcut = 3.57 cells. 

The simulations are started with the drop 0.1 R above the pool, where R = D/2 is the 
drop radius. Non-dimensional time is defined as = t/r, where r = D/V. The origin of 
time is taken when the undisturbed sphere would first contact the pool. The drop is kept 



Drop impact entrapment of bubble rings 



5 




Figure 3. Drop aspect ratio a ■ 



Dy/Dh vs. if. Th e bla ck oscillating line corresponds to the 
fitted equation (2.3). 



at a constant effective diameter D = 4.6 mm. Air has a viscosity of /i^ = 1.81 x 10~^ cP 
and density pa = 1.21 kg/m^. The liquid is water for both the drop and the pool, with 
viscosity fi = 1 cP, density p = 1000 kg/m^ and surface tension a = 72 mN/m. Gravity 
is included as ^ = 9.81 m/s^. We do not take into account the different properties of the 
bottom liquid in the simulation, and therefore do not include any Marangoni or variable- 
density effects between the two liquids. More details about the adaptive grid refinement 



can be found in Popinet (2003) and Thoraval et al (2012). 



3. Results and discussion 

3.1. Isolated bubbles and multiple bubble rings for a water pool 

We start by looking at the impact of the water drop onto a water layer. Figure |4] shows 
the early evolution of the outer neck contact of the drop with the pool. The contracting 
inner air disk is visible on the left side of the images in panels (d) and (e). Note that we 
are only looking at the early contact when the neck has not reached the size of the drop, 
as shown in the sketch in Fig. [ijb) . The radius of the neck in the last panel of Fig. ^d) 
has only reached 37% of the drop radius. 

Figure [4] shows that even for low impact velocities the neck region between the drop 
and the pool does not remain smooth and axisymmetric, but develops azimuthal undu- 
lations. For the lowest impact velocities these undulations have long wavelengths and 
do not entrap bubbles, see Fig. Qa-c). However, with increased impact velocity V the 
wavelength reduces and their amplitude grows more rapidly. In Fig. Qd) these undula- 
tions appear first in the second panel and grow in amplitude during the radial motion, 
but individual bumps saturate and are often being pulled back by surface tension. The 
shapes are irregular, but we can glean a characteristic wavelength from the third panel in 
Fig.|4]^d), giving A ~ 53 /im, corresponding to 73 undulations around the periphery. These 
undulations appear when the ejecta emerges, pulling local sheets of air under the ejecta 
on both or alternating sides of it. These local sheets can be pulled along with the ejecta 
base, with only occasional bubbles entrained, when these small azimuthal air discs make 
contact across the thin air layer, as is shown in a longer sequence of frames in Fig. |5j 
Individual bubble entrapments can also occur in the troughs between the undulations. 

In Fig. [6] we estimate the growth of the maximum undulation amplitude, measured 
between the troughs and peaks. The growth slows down with radial distance. For reference 
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Figure 4. Early contact of a water drop impacting a ^ = 250 /xm deep layer of water, (a-c) No 
bubble entrapment is observed at low impact velocities. An azimuthal instability develops on 
the ejecta, with a wavelength decreasing with increasing Re. (a) Re = 3610, We = 39, K = 2 360, 
a = 0.94, (b) Re = 4400, We = 58, K = 3860, a = 1.17, (c) Re = 8980, We = 90, K = 6640, 
a = 0.86. (d) For an intermediate impact velocity, individual micro-bubbles can be entrapped. 
Frames are shown at 1, 13, 18, 25 & 46 /xs after the first contact. Re = 11400, We = 394, 
K = 42100, a = 1.05. (e) For slightly higher impact velocity, the drop entraps one bubble ring 
and isolated bubbles, shown at t = 3, 11, 13, 21 & 40 /is. Re = 13300, We = 535, K = 61800, a 
= 0.98. (f) Multiple bubble rings. Frames are shown at about 5, 9, 13, 17, 21 and 32 /is after 
first contact. Re = 12 900, We = 506, K = 57600, a = 0.80. The scale bars are all 200 /im long. 
See also supplemental videos. 
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Figure 5. Details of individual micro-bubbles entrapments in the same conditions as in Fig.[4]^d). 
Frames are shown 2 /is apart. The top row shows the dynamics leading to the entrapment of a 
^ 10 fim diameter bubble. On the second row, two smaller bubbles of diameter ^ 4 and 6 /xm 
are separating from the right part of the edge. While the first one stays behind the edge, and 
can be seen in the last frame, the second one is re-absorbed into the neck in the last row. The 
scale bar is 200 /xm long. 



we plot a viscous length-scale 

L, = C^iy{t-to) (3.1) 

suggesting inertia and viscous forces both play a role. Here to is the time of first observed 
undulations on the front. 

Figure [5] shows that the characteristic azimuthal size of the undulations also grows 
during the radial motion of the front, but this is more difficult to quantify. 

It is curious that some air entrapment in breaking gravity waves has superficially 
similar appearance (Kiger & Duncan ( 2012| , their Fig. 11), but is clearly driven by a 
different mechanism and is three orders of magnitude larger in size. 

For slightly higher impact velocity, the entrapment of rather irregular bubble arcs be- 
gins. Figure [4]^f) shows up to 10 such partial rings. The average radial spacing of the 
adjacent bubble rings is ^ 26 /im. The air cylinders then break up into a row of bubbles 
through surface-tension driven Rayleigh instability. Bubbles are often shifted sideways in 
the azimuthal direction during the radial spreading (arrows in Fig. [4]^f) and supplemen- 
tary videos). For a stationary hollow cylinder of diameter djj in an inviscid liquid, the 
most unstable wavelength is = irdb/OASA. The characteristic time s cale of the expo- 
nential growth ~ exp(t/r(j) of the breakup is given by r^r = 1.22y^ pr^ja ( Chandrasekhar 
1961). The radii of the bubble arcs for a water layer, in Fig. [4]^f), are ~ 3 //m and they 
break up in about ~ 3 /is, which is 4rcr. 

Therefore, the first bubble-rings entrapment for water occurs around Re ^ 12 000 and 
K ~ 50, 000. These values are consistent with numerical results of Thoraval et al (2012), 
where no bubble ring entrapment is observed for Re = 10000 and K = 30000, and a row 
of bubble rings observed for Re = 14 500 and K = 74400. In the former case, the ejecta 
sheet is thicker, because of the stronger surface tension effects on the ejecta sheet owing 
to the lower value of the splashing parameter K. It is re-absorbed on the drop or the 
pool during the oscillations, and no bubble ring entrapment is predicted. However, in the 
latter case, at higher the ejecta sheet is thinner, and the oscillations entrap a row of 
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Figure 6. The maximum amplitude of the undulations of the front in the neck region, for water 
drop onto water pool. For th e sam e conditions as Figs.[4]^d) and [S] Data from three realizations. 
The solid line shows formula 13. 11 with C = 14. 



bubble rings at the core of vortex rings when it impacts onto and connects with the drop 
or the pool. 

However, the comparison of Fig. |4je) and (f) shows that more bubble rings can be 
observed at a slightly lower Re and K. This suggests that the Re number of the impact 



is not enough to characterize the bubble-rings entrapment. We will show in ^3.4 the 
critical effect of the drop shape. Moreover, the azimuthal instabilities also affect the air 
entrapment, and individual bubble entrap ments have been observed at slightly lower Re 

in Figs, ^d) & [S] 

Note that in the work of |Castrej6n-Pita et al. 1( 2012), the vortex street is observed for 
conditions similar to Fig. 4(f) of [Thoraval et al. (2012), where no bubble-ring entrapment 



was predicted. Considering the large diameter of the drop they are using, it is not clear 
if the bubbles they observe are part of a bubble ring or isolated bubbles. They are also 
looking at a larger view and longer time evolution, that could be out of the field of view 
used in the current investigation (see perspective in Fig. [T]). However, their alternating 
vortices are a new observation, clearly different from the isolated vortex rings produced 
by much lower impact velocities, see [Peck fc Sigurdson | ( p"994[ ) . 

3.2. Bubble rings for miscible liquids 

The bubble ring entrapment becomes more regular for the impacts onto ethanol and 
methanol pools (Fig. [t]), perhaps due to the lower surface tension of these liquids (Ta- 
ble [T]). Contrary to the impacts onto water films, no azimuthal instability develop on the 
neck of the ejecta, which remains perfectly smooth in Fig. [7|^a-c). At higher impact ve- 
locities bubble rings are entrapped. Figure [7|(d) shows at least 10 bubble-rings. Many of 
them are entrapped axisymmetrically over the entire image view, which can span around 
90^ angular sector. 

In some instances thin ribbons of air are entrapped and subsequently breakup into 
sub-rings and thereafter bubble rings, as highlighted by arrow in Fig. [7|^d). The frames 
in Fig. [7|^f) detail this sequence of air entrapments (green arrows). The air sheet first 
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Figure 7. Early contact of a water drop impacting a pool of a miscible liquid (ethanol or 
methanol), (a-c) No bubble entrapment is observed at low impact velocities on a thin film of 
ethanol (6 250 /xm). (a) Re = 3610, a = 0.94, (b) Re = 4400, a = 1.17, (c) Re = 8980, a = 
0.86. (d) Bubble rings for a water drop impacting onto a methanol layer (Re = 12 900, a = 0.80, 
S = 50 /im). Frames are shown at 5, 9, 11, 15, 20 & 33 /is after first contact, (e) First oscillations 
of the ejecta sheet, followed by entrapment of bubble rings, for a film of ethanol (Re = 13 300, 
a = 0.98, 6 = 250 /xm). Azimuthal instabilities appear in the ejecta sheet, before its rim detaches 
in a liquid toroid. The first 8 frames are shown 4 /is apart, and then 20 /xs. (f) Bubble entrapment 
dynamics for impact on a methanol film (Re = 12 900, a = 0.80, ^ = 50 /xm). The first frames 
shows the entrapment of a superposition of air sheets, later breaking in patches and then in 
micro-bubbles (green arrows). Three bubbles are identified in the third frame by red, black and 
blue arrows. Bubble arcs with legs in the radial direction are identified by the yellow arrows. 
Frames are shown 3, 6, 15 and 60 fis after the first one. The scale bars are all 200 fim long. 
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Figure 8. Larger view of the drop impact onto ethanol, at times t ~ 50, 150, 250 & 1150 /is 
after first contact {Re = 14 500, a = 0.94, 6 ^ 125 /xm). Frames taken from videos using Photron 
SA-5 at 10,000 fps, with a 1 /is exposure time. The central bubble has drifted out of the image 
in the last frame. The scale bar is 1 mm. 



breaks into smaller patches, later contracting into bubbles. In the third frame, some 
bubbles are observed at the same radial location as another air cylinder. This supports 
the mechanism that air can be entrapped both above and below the eject a sheet, as was 



shown in the numerical simulations of Thoraval et al. (2012). 

In a similar way as for water pools, Fig.[7][f) also demonstrates strong sideways motions 
of bubbles. Three bubbles are identified in the third frame by red, black and blue arrows. 
Their corresponding location is marked by the same colored arrows in the last frame. 
Comparison between the black and red arrows shows that this sideways motion can be 
of different strength for adjacent rings. This bubbles motion in the azimuthal directions 
results in their clustering at isolated locations. Bubble arcs with legs extending in the 
radial direction towards the neck, are identified by the yellow arrows. 

Figure |8] shows a wider view of an impact taken at a lower frame rate of 10,000 fps but 
larger pixel area of 896 x 848 px, using the Photron SA-5 CMOS camera. Each frame is 
frozen with a 1 /is exposure. The 100 /is interframe time only shows us snapshots of the 
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Figure 9. (a) First onset of bubble rings for water drop impacting onto a water, ethanol or 
methanol layer of different depths (see also Fig. 13 'b)). (b) Spacing of adjacent rings Ai?, vs 



radial distance R from the impact center. Data from 8 different realizations for a methanol pool 
{Re = 13 300, a = 0.98, S = 250 /xm). 



phenomenon, putting the earlier figures in perspective, with most of the earlier sequences 
occurring before the first image. In the second frame multitude of splashed droplets 
appear from underneath the shadow of the drop, with some larger droplets planing on 
the pool surface, leaving behind narrow capillary wedges. The first panel shows a smooth 
central regions, followed by a convoluted interface, suggesting the stirring by the three- 



dimensional vortical structures (see ^3.7). Similar stirring can be inferred from the side 
shadowgraph imaging in Castrejon-Pita et al. (2012) (their Fig. 4). The final panel in 



Fig. |8] shows numerous isolated bubbles which have been redistributed by the vortical 
motions. The bubbles are mostly concentrated within the mushroom-like remnants of the 
vortical structures. 



3.3. First onset, number and spacing of bubble rings 

The first contact entraps a central air disk and forms a rapidly expanding outer liquid 
edge. Bubble rings are then formed as observed above. Figure [9]^ a) looks at the radial 
location where the first bubble ring is entrapped. We normalize the radius of the first 
ring Ri with the horizontal drop radius when it first contacts the pool, Rh = Dh/2. The 
data shows large spread, but an overall trend is for the onset to occur earlier for larger 
impact Re. The lowest entrapment radius is 0.18, similar to the 0.2 limit observed by 



Thoroddsen 



( 2002| ) for the onset of the ejecta sheet. This onset radius of the ejecta is 

also in agreement with the inviscid numerics of "Weiss & Yarin (1999). 

Figure |9jb) shows the distance between the adjacent bubble rings measured for numer- 
ous identical impact conditions. The spacing of the rings tends to increase with distance 
and the entrapped bubbles become larger. 



3.4. Effect of pool depth and drop shape 
To ascertain the influence of the pool depth we systematically vary the layer thickness ^, 



from about 25 /im to 1 mm. Figure lO'a-c) compares the bubble rings for the 3 smallest 
pool depths S. The ring structures are qualitatively similar in all cases, but the shallowest 
pool shows the earliest and finest bubble rings, some of which are sub-pixel in diameter. 
The second ring in Fig. 10 ^a) allows us to measure the separation of micro-bubbles, giving 
A = 8.8 /im, suggesting a diameter of the original air torus dtor - 1-4 /im. The earliest 
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Figure 10. Effect of pool depth on the bubble rings for water drop impacting onto methanol, 
with S ^25 (a), 50 (b) and 75 /xm (c) {Re = 12 900, a = 0.80). The frames are all shown 24 /xs 
after the first contact in (a-c). (d) Deeper methanol pool with 6 ^ 500 /xm. Frames shown at 10, 
14 & 28 /iS after first contact. The scale bars are 200 /xm long. See also supplemental videos. 



ring appears even smaller, arrow in Fig. 10 ^a). Figure [TQ[d) shows fewer by qualitatively 
similar ring entrapment, for much thicker layer. 

Figure [TT] shows numerical results for three different pool depths, where the two-liquid 
interface is highlighted by coloring the drop and pool differently. It shows the drop 
penetrating further into the pool for larger S. The vortex street is therefore constrained 
in a shallower region for shallower pools, and therefore develops more horizontally. This 
constraining effect increases the maximum liquid velocity by as much as 20%, as shown 
in Fig. 13 ^a). The first bubble ring is entrapped at Ri/R of respectively 0.37, 0.31 and 
0.28 for pool depths S = 800, 200 and 100 /im respectively. This confirms the previous 
experimental observation of Fig. 10 that earlier rings are observed for shallower pools. 
Figure 13 'b) also shows this experimentally in a more systematic way, but the difference 
is not very pronounced. The radial location of this first entrapment is also consistent 
between experimental and numerical observations, even though numerical simulations 
only considers one liquid. 

In all three cases, we can see in the numerics that the central air disk punctures at the 
center during its contraction into a central bubble, thus forming an air torus. This is also 
observed in some of the experiments, as seen in Fig. [T2j where an air torus is formed, 
which later contracts into one bubble. 

We observe experimentally that the most robust bubble rings are produced by a fiat- 
bottom drop in Fig. [4]^e) & (f) and Fig. [7|^d-f) with fewest bubble rings in 7(e), which 
is more spherical. The largest number of rings is also produced by such oblate drops 
(Fig.|4jf) and Fig.[2TJb)). This is consistent with the numerical results of Fig. 14, showing 
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Figure 11. Effect of pool depth on the vortex street at t* = 0.0526, for Re = 13 700. Pool depth 
S = 800 /xm (a), 200 /am (b), 100 /am (c). The top of the 3 images are at the same location, 
relative to the original pool surface. The pool depth is larger than shown in the image in (a), 
but is completely included for (b) & (c), where the bottom is indicated in gray. The maximum 
and minimum level of refinement in the domain are respectively 12 500 and 778 cells per drop 
diameter and Dcut = 1.3 /xm. The bar is 200 /xm long. We can observe that one part of the ejecta 
sheet is climbing on the drop in (b), while the ma in ejecta sheet continu es to emerge below. 
Th is is similar to what was observed numerically bylThoraval et al. (2012) and experimentally 
by Zhang et al. (2012). The formation of this higher part of the ejecta sheet can be observed 
in the Supplemental Video. It then merges with the drop in this case. 



Figure 12. Contraction of the central air disk for a water drop impacting on a thin water film 
{Re = 11600, a = 0.86, S = 250 /xm). The convergence of capillary waves punctures the center of 
the disk, producing an air toroid in the third frame of the second row, later contracting into a 
spherical bubble. Similar formation of a central air toroid can be observed in the supplemental 
video of Fig. p[e), as well as in other water drop impacts on ethanol or methanol films. Frames 
are shown 7.1 fis apart, with an exposure time of 1 fis. 
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Figure 13. (a) Maximum velocity Umax in the liquid from the numerical simulations in Fig.|ll[ 
for 3 different pool depths, (b) Experimental observation of the first bubble ring entrapment 
radius (Re = 12 900, a = 0.80) for different pool depths. 




Figure 14. Effect of drop shape on the early dynamics at t* = 0.0406, for Re = 12 900 and 
5 = 800 /xm. The drop is modeled as an ellipsoid of revolution, where a is (a) 1.29, (b) 1, (c) 
0.79 . (a) and (c) correspond to the maximum horizontal deformations of the fitting equation 
(2.3), with a = 0.162, keeping the same effective diameter. The maximum and minimum level 
of refinement in the domain are respectively 5240 and 655 cells per drop diameter, and Dcut = 
3.1 /xm. The scale bar is 200 /xm long. 



a larger number of rings for the oblate drop. It even suggests that a prolate drop could 
completely suppress the bubble ring entrapments for the same effective diameter, as is 
shown in Fig. [l4]^a). 

3.5. Edge breakup and splashing 

Numerical simulations have shown that the ejecta sheet can impact alternatively on the 
drop and the pool during the vortex shedding. The tip of the ejecta can thus detach into 

) . Such tori are 

highly unstable to Rayleigh instability and break rapidly into splashed micro-droplets. 



a liquid torus exiting the neck region at high speed ( Thoraval et al. [2012 
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Figure 15. Details of the ejecta sheet breakup; for a water drop impacting a thin film of ethanol 
(Re = 11700, a = 0.79, ^ = 75 /xm). The ejecta sheet starts to puncture on the second frame of 
the second row. The growth of the holes leave tendrils connecting the neck to the liquid rim. 
The scale bar is 200 //m long. Frames are shown 2 /xs apart. 



However, this earliest splashing of micro-droplets by axisymmetric breakup of the ejecta 
sheet had not been observed previously in experiments. 

By looking carefully at Figs. [7]^d-f) & [sj we can identify this early splashing by the 
breakup of the tip of the ejecta sheet after the entrapment of a few bubble rings, as was 
suggested by the numerical simulations. The liquid toroid in Fig.[7|^e) separates at ^ 25 
/is after the first contact, and a velocity of 20.5 m/s, which corresponds to 7.1 times the 
impact velocity. The tip velocity of the ejecta sheet in the numerical simulations at the 
same non-dimensional time, is 6.2 for Re = 12 900 (Fig. 14 ^b)), and 7.4 for Re = 13 700 
(Fig. [iTJb)), which is in excellent agreement with the experimental observations. 

Figure [15] shows more detail of the breakup of the edge at a slightly lower Re. The ejecta 
sheet breaks via holes puncturing behind the rim. The thicker ejecta rim is therefore left 
connected to the neck by liquid tendrils, but subsequently becomes fully detached. 

To remove the ambiguity of the bottom view, we have also looked at this early splash- 
ing from the side above the free surface. Figure p^a) confirms the bottom view images, 
showing the ejecta sheet emerging from the neck, puncturing behind the rim and sepa- 
rating a liquid toroid from the neck. Liquid tendril are also observed in the fourth frame, 
and are slingshot ahead of the rim, creating the protrusions observed in the last frame, 
similar to the last frame of the first row of Fig. ^e). The slingshot of the broken ejecta 
sheet can also be observed in numerical simulation, as in supplemental video of Fig. fl4ja) , 



and is similar to the slingshot mechanism described in [Thoroddsen et al. (2011). 



Two consecutive liquid rings are observed in Fig. 16 'b). It also shows the emergence of 
a greatly disturbed ejecta sheet after this early splashing. A larger bottom view confirms 
this side view observations in Fig.[T6|c). This mechanism, of a detachment of a thin torus 
of liquid, explains the synchronized emergence of uniform sized micro-droplets observed 
ahead of the main irregular ejecta sheet, see Fig.[T6jd) as well as Fig. 2(c) in Thoroddsen 
(2002). The irregular sheet is clearly shown in Fig. 8(b). 

The splashing of several liquid tori is consistent with numerical simulations showing 
that the ejecta sheet can breakup during the successive impacts on the drop and the pool. 
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Figure 16. Detachment of liquid toroids by breakup of the eject a sheet. Scale bars are 500 /xm 
for (a,b) and 1 mm for (c,d). All images have an exposure time of 1 /is. (a) Side view of the 
liquid toroid detachment for a thin film of ethanol (Re = 12 600, a = 1.11, S = 125 /xm). The 
eject a sheet punctures in the third frame, breaking into tendrils. Those tendrils are slingshot in 
front of the liquid rim which has not yet broken, as observed in the last frame. The frames are 
shown 10 /iS apart, (b) Side view of two liquid toroid detachments for a thin film of methanol 
(Re = 12 600, a = 1.11, 6 = 125 /xm). The frames are shown 10 /xs apart for the first 3 frames, 
and then 80 /xs. The splashing of liquid toroids is followed by the emergence of an irregular 
ejecta sheet, as can be seen also on the larger bottom view of Fig. [s] (c) Larger bottom view 
of two consecutive rings of liquid droplets and toroidal sections detaching for an ethanol pool 
(Re = 11600, a = 0.86, 6 = 125 /xm). (d) First droplets emerging from underneath the drop in 
the same conditions as Fig. [S] Image difference between two adjacent frames, to highlight the 
splashing droplets. Shown in inverted gray-scale. The exposure is 1 /is long. 



Supplemental video of Fig. [T4[a) shows such a case where the ejecta sheet breaks first 
by climbing on the drop and then impacting on the pool, thus creating two consecutive 
liquid tori. 

After the first ring entrapment, regular spanwise instabilities can appear in the ejecta 
sheet, as is clearly seen in the second panel of Fig.[7|^d), as well as Fig.[7|^e) & (f). The fine 
azimuthal breakup when the ejecta bends and impacts onto a pool have been reported by 
Thoroddsen et al. ( 2011[ ) (their Fig. 5) and may be of similar origin. Furthermore, the 



early appearance of similar azimuthal instabilities have also been observed by |Thoroddsen 



et 



( 2012, ) in a free-surface cusp, which is formed during a drop impacting onto a solid 
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Figure 17. Entrapment of a large air cylinder for a water drop impacting a thin film of methanol 
{Re = 11400, a = 1.05, ^ = 50 /xm). The small bubble at the bottom right of the long cylinder 
make a complete rotation around the cylinder between frames 4 and 12. The frames are shown 
4 /is apart. The scale bar is 200 /xm long. 



surface. Numerical simulations show that the ejecta sheet breaks when it impacts on the 
drop or the pool, by stretching between the new connection and the faster rim. This 
instability is therefore consistent with the impact of the ejecta sheet on the drop or the 
pool. Similar breakup of a liquid sheet by stretching was also observed experimentally 



by Roisman et al. (2007) for spray impacts. 



3.6. Vortex shedding and rotation around bubble rings 

The difference in refractive index between the drop and the pool (see Table [T]) allows us 
to visualize vorticity structures inside the liquid. As the coherent vortices bend and wrap 
up the interface between the two liquids, a dark line can be observed at their edges with 
our back-light imaging setup. 

Numerical simulations have shown that the first oscillations of the base of the ejecta 
sheet have a smaller amplitude and do not entrap any bubble rings, see |Thoraval et al~| 
(2012) (their Fig. 4(c,f,g)) and our Figs. 11 & 14 This is consistent with our experimental 
observations of Figs. [Tj^e), 15, 19]^a) & 21, where dark arcs form before the first bubble 
rings. They show the shedding of vortices from the neck before the start of the bubble-ring 
entrapment. 

As the neck moves outwards radially, the angle between the pool and the drop becomes 
larger, and bubbles rings are entrapped, as observed above. Numerical simulations have 
shown that these rings are often entrapped in an alternating way at the top and the 
bottom of the ejecta sheet. At the same time, they shed vortices of alternating sign in 
the liquid. Bubble rings are therefore entrapped at the core of vortex rings. Dark lines are 
indeed observed experimentally around the bubbles, supporting this shedding of vorticity 
scenario. The rotation is also made apparent by the dynamics of a micro-bubble rotating 
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Figure 18. Tracking of a 2 /im particle rotating around a vortex core for a drop impacting on 
a pool of ethanol {Re = 11 400, a = 1.05, 6 = 125 /xm). AL is the distance to the initial positions, 
minus the mean radial translation component identified by a linear regression. The period of 
rotation is here about 18 /xs. 



around a larger bubble cylinder in Fig. 17 We can obtain an estimate of the rotation 
speed in the vortex rings by tracking a 2 jam particle seeded into the pool liquid (Fig.[l8|, 
giving a rotation period of 18 /is in this example. 

The dynamics of the air entrapped in this vortex street is also affected by the rotation 
around it. One can expect the rotation to delay the capillary breakup of the air cylinders 
( Rosenthal 1962; A shmore fc Stone ||2004] [Eggers fc Villermaux ||2008[ ). This is indeed 
evident in Fig. [T7| where the cylinder of diameter about 29 fim can be observed for 
more than 94 /is before breaking, corresponding to t/ra- - 7 based on the methanol 
properties. The air cylinder also elongates by about 25% between the first and the last 
frame, consistent with the theory that it resides inside a vortex. The relative motion 
of the bottom tip of the air cylinder and the closest micro-bubble below shows that 
the stretching due to the radial motion cannot account for this elongation. Even smaller 
cylinders can be stabilized, as observed in Fig. 10 'c). Two air tori are formed next to each 
other, with similar diameters dtor - 8 /im (see arrows in the figure), but break up at very 
different times from formation, of t/rcr - 4 and 12, based on the liquid properties of the 
methanol. These delayed breakups show the strong stabilization effect of the circulation 
around these air tubes. The breakup wavelength is also larger, as the theory suggests 
( Ashmore fc Stone ||2QQ4l ). 

The row of vortices shed in the liquid can also interact with adjacent ones. In some 
realizations two closely entrapped air tori rotate around each other, with the line of small 
bubbles rotating around the bigger one. Figure 19 ^b) shows such a sequence, where we 
track one rotation, which takes T = 104 /is. Numerical simulations show that vortices of 



different strength are created at the top and bottom of the eject a sheet ( Thoraval et al 
] |2012 ). The longer time evolution shows that this difference can make the bubble rings 



rotate around each other while translating (Fig. [T9[c) and supplementary videos). These 
dynamics are consistent with the experimental observation of the rotation around bubble 
rings. 

Some of the bubble tracks simply translate past each other during their radial mo- 
tions, for example visible the videos accompanying Figs. 19(a) fc 21(a). The bubbles 
are initially sitting at slightly different depths, and this difference in vertical location 
is amplified by the vorticity, as shown above. Moreover, as the bubble tori break into 
bubbles, their vertical width will slightly increase, sampling larger mean shear. These 
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Figure 19. (a) Bubble rings formed during the impact of a water drop onto a shallow pool of 
ethanol (Re = 13 300, a = 0.98, 6 = 250 /xm). The frames are shown at 8, 20, 38, 66 & 112 /xs after 
the first one, showing a total of 7 bubble rings, (b) Close-up view of a line of fine bubbles which 
circulate around another line of slightly larger bubbles, from the video in (a). Total duration of 
these frames is 104 /is. See also supplemental videos, (c) Numerical simulations of a water drop 
impacting on a thin film of the same liquid (Re = 13 800, a = 1, 6 = 800 /xm). This focused view 
of the interface at t* = 0.481 shows the rotation of a pair of bubble tori at the core of adjacent 
vortex rings. The scale bars are all 100 fim long. 



relative translations of bubbles could therefore result from a vertical mean gradient of 
horizontal velocity within the pool depth. 

3.7. Three-dimensional instabilities 

The axisymmetry of the impact is rapidly broken by different instabilities (see Figs. |4] 
&[7|). For a water pool, we have observed undulations develop in the neck as soon as 
an ejecta forms (Fig. |4|. This leads to entrapment of isolated bubbles (Fig. [4] & [s]) and 
creates less regular bubble rings compared to ethanol or methanol pools. Even in those 
lower surface tension liquids, the most regular rings appear at smaller entrapment radii, 
where the ejecta sheet has not broken yet and remains axisymmetric. 



We have already suggested in ^3.5 that some of the azimuthal instabilities come from 
the breakup of the ejecta sheet, when it bends and touches the drop or pool surfaces. Over- 
turning gravity waves also rebound and destabilize underlying vortices, but at a much 



larger scales than herein, see|Watanabe et al. (2005). Figure 15 clearly demonstrates the 



effect of isolated neck disturbances on the bubble entrapment. A small perturbation is 
visible in the first frame and develops in time. The two first bubble rings which form in 
the following frames are broken at the same azimuthal location and the ejecta ruptures 
there first. 

Three-dimensional instabilities also develop inside the liquid, and are made apparent 
by the difference in refractive index. Radial lines are visible in Figs. [T^e), § ^^), [2^ 



& 21 They show the formation of streamwise vortices between the primary spanwise 



vortex rings, which often reach to the free surface in the neck. For the lower Re cases 



(Figs. 20), isolated streamwise vortices are observed. In Fig. [2Q[a), they appear in pairs, 
at the same location as a front perturbation. The lower one in the image starts at an 
isolated location on the side of a vortex ring, where it entraps one micro-bubble. Two 
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Figure 20. (a) Roll- up of isolated streamwise vortex pairs for a water drop impacting a film of 
ethanol seeded with 2 /xm particles (Re = 11400, a = 1.05, S = 125 /xm). The first images show 
that the vortex pairs are starting at the same azimuthal location as local disturbances in the 
front. The bottom vortex pair first entraps a bubble, later splitting in two. Frames are shown 
3, 6, 9, 12, 17, 23, 27, 35 /xs after the first one. (b) Roll-up and sideways motion of streamwise 
vortices for a water drop impacting a film of ethanol seeded with 2 /xm particles (Re = 11400, 
a = 1.05, 6 = 125 /urn). Frames are shown 3, 10, 18, 32 & 47 /xs after the first one. The scale bars 
are 200 /xm long. 

lines are then visible on each side of this initial entrapment and extend up to the front. 
This suggests that the vortex pairs arises from the same vortex loop, rolling- up around 
the vortex ring. The connection between the two streamwise vortex lines should then 
form a vortex loop with one section in the azimuthal direction, near its origin around 
the bubble. The presence of a strong vorticity around this bubble is demonstrated by 
its breakup into two smaller bubbles (in the 4th panel of Fig. [2Q|^a)). Moreover, the 
later dynamics shows the secondary streamwise vortex tubes roll up around the primary 
spanwise vortex ring. Figure |20]^b) shows a similar case where streamwise vortices roll- up 
around a vortex ring. This roll-up can be identified by following two micro-bubbles at 
the core of the vortices (black arrows). 

Interesting parallels can be made with similar three-dimensional instabilities occurring 




Figure 21. Bubble rings and vortex entanglement {Re = 12 900, a = 0.80). (a) Frames shown 
are 8, 22, 41 & 58 fis from first contact on ethanol (6 = 250 fim). (b) Careful examination of the 
fine bubbles, show 14 separate bubble rings, which are pointed out by the arrows, (c) Closeup 
of the vortex tangles at t = 58 /xs. The scale bars are 100 /xm long. See also supplementary video. 



in the wake of a cylinder (Williamson 1996). The local Reynolds number at the base 



of the eject a sheet Rcb will be affected by both the radial velocity of the neck and the 
velocity within the liquid. Rei, can be increased both by increasing the impact velocity or 
by using a more oblate drop. Indeed, as the impact Re increases, the eject a sheet emerges 
from a faster moving neck ( Josserand fc Zaleski ||2003 Thoraval et al. |[2012 ), and at a 



higher velocity (Thoroddsen 2002 Josserand & Zaleski 2003), thus leading to higher 



Rc}). A flat bottom drop also geometrically generates a faster moving neck, and produces 
larger velocities in the liquid, see FigJlSja). Both effects lead to a larger concentration 
of streamwise vortices (Figs.[7|^e),[8l|19[a),[2Q|fc[2T]). This is similar to the onset of three- 



dimensional instabilities of the vortex street behind a circular cylinder (Williamson 1988 
|199 6), where finer streamwise vortices and a smaller spanwise wavelength are observed 
at higher Re. Similar vortex loops are also observed in both cases, as described above. 
Moreover, rapid motion of the bubbles in the spanwise direction along the vortices is 
observed in our experiments, see Figsjljf) and [7|^f). Similar lateral motion was also 
observed behind "vortex dislocations" in the wake of a cylinder ( Williamson [1992 , 1996| ). 



4. Conclusions 

Observations from below the impacting drop have herein demonstrated that the mech- 
anism suggested in Thoraval et al. (2012) does indeed entrap micron sized air tori. The 
oscillations of the eject a sheet can thereby entrap a row of bubble rings. The vorticity 
entering the liquid during those oscillations and bubble entrapments is then destabi- 
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lized into complex 3D structures. The combination of azimuthal instabilities and vertical 
oscillations produces a large range of new bubble entrapment scenarios. 

Besides imaging the formation and breakup of partial bubble rings, of equal significance 
is our observation that the outer neck is unstable to azimuthal undulations, for water on 
water impacts at even very moderate Reynolds numbers (Figs. 4(d) & 5). This poses a 
challenge to theoretical and numerical studies, which invariably assume axisymmetry. 

We note that the bubble rings observed herein differ in fundamental ways from the 
Oguz-Prosperetti rings ( Oguz fc Prosperetti |1989| , as the base of the ejecta is not driven 
by surface tension, but rather by the impact pressure. This high localized pressure is 
indeed the mechanism responsible for driving out the ejecta sheet. 

However, the details of the air entrapment and its dependence on the impact conditions 
is still not clear. Large scope exists for further work, as the present study perhaps raises 
as many questions as it answers. What role do Marangoni or Rayleigh-Taylor play in the 
two-liquid dynamics? The large parameter space of other liquids needs to be studied. 
Even for the same liquids, in the drop and pool, the interplay between Re, We and 
a which allows entrapment, or preserves an extended axisymmetric ejecta, remains to 
be determined. The three-dimensional instabilities of the vortex street also need to be 
studied in more detail and compared to the well-known instabilities of the cylinder wake 
and the shear layer. However, in the vortex street observed here, the vortices are shed 
from a deformable free-surface, adding to the complexity of the analysis. The influence 
of the shed vorticity on the dynamics of the neck, both vertically and in the azimuthal 
direction, should also be added in the stability analysis of splashing. 
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